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plus ix 

The iSilrlwind I ooopiiter, liWX, at KIT le a lilgh*8peed eleo- 
tronio digital conpater* It will be used in a wide variety of applloatione 
Involring coapat&tlonB f or aatbeaiatleal and engineeriag problems, aooountingp 
Btatietioal analTfle, and slnolatlon and control prooeeeee* The prepare* 
tlon of eaeh problema for ffKI le essentially a setter of eiqpreesing tbe pro- 
cesses to be perfoiaed in terse of a list of tbe oosputer^e basic operations* 
Such a list Is called a progren, and there are a sailsua of 32 basic opera- 
tions, each of which WVI performs automatically upon receiving a single in-' 
struction from a program* WWI also automatically seqaeitbes the operations 
listed in a program so that they are performed one after the other* 

The terminology and concepts of prograoBlng for fWI and of IVH 
itself, essential to an understanding of programming for the computer, are 
presented here In the eeetion on Eseentiale of Progrmmmlng and in the Ap« 
pendiz* A progreeeively developed eet of esamples of programs Is ihduded 
to indicate hoe programa are developed, to illustrate the nature of various 
nays of handling dixferent prooeases, end to suggest the flexibility of pro- 
gramming techniques* Use is made of flow diagrams and of subroutine tech- 
niques, both in the aaalysit of problems and In the preparation of programse 

Once the prograMer understands the basic IWI operations and 
how programs are developed from them, the programming of a problem reduces 
to an analysis of the original problem* This analysis la the same, of 
course, as if the actual computations were to be made hy lass automatic 
operations* 



<^ 
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FOREWORD 

This report it an introduction to prograiunlng for tbe 1lhirl<» 
wind Conputer, WWI, at MIT* 

It begone with a general atatenent of what WKl it and what imi 
oan do* This atatenent ie intended to suggest both the possible areas of 
applioation and the linitations of a high-speed digital oonputer suoh as 
wni* 

Fart I presents essentials of WWI progransingo 

Part II presents exaoples of WKI progranso The exanplee are 
ohosen to illustre^ > bow progress are developed and to suggest various 
programing techniqueso 

Part III is an appendix which includes a discussion of nuii» 
bars in wm and the Ml operation code* 
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INTRODUCTION 

URAT WWI IS 
HEAT mi COES 
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WPAT WH? |§ 

nwi 10 A hlgfa«gpeed fl^g^Wtrf? 4i£i^ SSi&H$AX* I^ is ^ nP^^ 
and verMtilo tool for prooeesing lofonBAtlon* Built into it «ro olootronie 
eirouitfi for aatomtioally Mqvwnoing and pox^foraiing Mveral baalo opara'^ 
tionft* Coflrplioatad prooateas nay ba da^alopad froa saquanoaa of thaaa ba- 
ale oparationeo j^ropriata inatraotione for tba parformanea of ttaata opar- 
ationa ara praparad in tha fom of progranfi* Prograaa ara raad into and 
atorad in imi intamal atoraga* WWI than oarriaa out thaaa instruotiona* 

Tha high apaad of oparation of IVI ia aada poasibla bgr tha ii*9 
of alaetronio eirauita inataad of aoving partt uead in alotar davioaao Op~ 
aration tina ia aaaaarad in Billlontha of a aaoond* An avaraga of tan thou- 
aand oparationi aueh aa addition, oultiplieation, and division ara parfonad 
in ona aaoond* At thia rata, in 15 ainutat WKI parfonna arror*frea calcula- 
tions which would raquira 15 yaara of oontinuoue hand oaloulationo 

All Wfl^'e circuits work on tha basis of thair ability to parfom 
tha aii^la taak of distinguishing bataaan tha prasanoa and^ abaanoa of currant 
flow* Tha digits 1 and ara assigned to tha two atatass currant on and 
currant off* Tha oonputar doas all its oaloulating bgr counting aith thasa 
tao digita alona, in tha binary systan, just as wa ara accustomed to calcu** 
lating vith tan digits in tha daoiaal syatev* 

Since KV>I fiSBSU «i^ discrete digits, it is called a digital 
cosqputar* Thia diatlnguiahas it froa analoif devices vhicb aeasur^ eontlnu** 
Ottsly changing quantities and give analogous physical, rather than mmerical, 
indioationa of aagnitudeo 
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Upon receiving appropriate Inttruotiont fron e previouely pre* 
pared progran WI atttoeatioally perfomt fequenoet of its baaic oparatione 
tery rapidly p ao tbats 

1) It oan eave nany hours of labor on routine cenputational 
tatka* 

2) It Bakes praetioal the perfomence of many BatheBatioal and 
0tati8tioal manipulationt which would be otherwiaa too 
oofitly or tine coneuaiing* 

3) It nakee poaaible the autonatic oontrol of eany proceatee* 

The following drawing auggeste a few of the poeeible applicationa 
of high-Bpeed digital coiiputera* 
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Ab with all tools 9 tho offootlTO utage of WICI z^ulrot a ItilX vp^ 
proeiatlon of what it it dotlgntd to lie* Baaioally it oan in eiagla opora- 
tionat 

1) Stora and timnafar infomation within itaolf* 

2) Parfom tha arithnatio oparationt of additiont snbtraetion, 
Bultiplioationy and diwiaiono 

3) Distinguish batwaan poeitiva and nagativa mabara and taka 
ona of two altamata ooursas in a prograa on tha haais of 
tha distinotion* 

A) Sand alaotrioal polsas to various axtamal davieas and con- 
trol thair action accordingly* 

It can automatically parfora thasa singla opavations in aaquancao 
Sequancas of thasa sinpla oparations nay ba conbinad by tha prograanar 
into vary langthy and conplicatad prooaeaas* 

Bacauaa thara is a considarabla sisilarity batwaan oparations in 
such aachlnas as Vil^T and thosa in a siapla narvous systan, and bacausa 
thasa oparations arc parfomad so rapidly, it has bacona popular to think 
of such aaohinas as giant brainso This analogy is not vaty halpflil to ona 
who is. preparing a program for tha ooaputarl His job is to spacify avary 
operation in tha proease which tha oomputar is to carry out* This naaas 
that ha must first oomplataly analyse tha process, using general terms in* 
stead of specific walues* Than ha must prepare a program in coded fonn 
(which describes the procedure step by step), to which the oeaputer is built 
to respondo Ihan tha oomputar is given this program it prooeeda to carry 
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it out frith cpeoifio miacrioal valuee, also supplied either hgr the pro* 
gremer or hgr extemel devieeso 

Thie relationehip of progrenRer to WKI detenree graphic il* 
laetretion* While etadjring the following drawing, inagine the operator 
of the defik oaloulator ae the prograaBer for Whirlwind I* Taagine hia to 
the left of the heavy arrow labelled "problem to be done"* In preparing 
the progrev of the problen, he nuet go through the wety e&ae analyeia whioh 
he goee through for the eame problea when using the desk oaleulatoro 

All oaloulatlcne handled b^ T.WI can also be handled on a desk 
oaleulatoro I'Wl^B usefulness lies in its speed in oarrTing out processes 
involving conputationo This nake's the perfoniiiance of aany processes eoon- 
onically Justifiable for the first tine; further it aakes possible auto^ 
■atic control of certain processes whioh, because of the speed required , 
could be handled in no other way« 

Of course progrfis preparation time nust be included In aqgr esti- 
■ate of MI*iJ cost In solving a particular probleno For this reason muaiy 
probleas whioh oaa be eolved nanually in a reasonable tiae, end which sseed 
not be solved often, aay not be wultable for solution on a large conputero 
On the other hand, Bany very siaple routine tasks can be econoaically 
handled by ^.I because they can be progranmed once and for al'*^ ' »«;a 
elioinating the need for separate analysis each tine the task ^s perforeedo 
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A* WHiar A fBOGSUM IS 

A progruD le • ••quwaot of InBtruotlona and nuMrleal ^raXaet, 
in coded font* It isstfuete tho conpator to porforn a partienlar prooeas, 
ono oporatlon at a tlsa* 

Bafora a prograu la parfomad^ it la atorad, in appxopriata 
form. In tha atoga^ tmit* It la than astacutad bsr tba oonbinad oparation 
of tha oontrol and arithaatle unlta* 

WIKI oan dlatlngnlah taatwaan and paxf oza 32 haalo kinda of op- 
amtlooa* Thaaa oparatlona ara daaoribad In "Tha Ordar Coda" whloh la In- 
olttdad In tha Appandlx (D)* 
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TlM iBrtnrtlW *i^ nytftf in a progTMi are both oallad 
words* IfiBdi are atorad la aaparata ragietara of tba atoraga tmlt* 
Each ragiatar la nnbarad, ani tba nwibar daalgaating a ragiatar ia 
oallad ita addraay * 

A prograMMd nihtX apaoif iaa tba alga and nagoituda of tba 
nunbar* 

An Inatnietion la a ooablnatlon of ona of tba 32 oparatlona 
plua tba addraaa of tba atoraga ragiatar ablob oontalna tba word to ba op* 
aratad on* 

Aa an axaapla of tba aamct aaaolng of an inatruetlon^ auppoaa 
tbat at aoaa tlaa In a prograa It la aaoaaaary to aubtraet -** 1/2* Tba 
progrunar aould Inelnda two words In tba prograa for tbla purpoaa* Ona 
would ba tba nuabar * 1/2* Tba otbar would ba tba Inatxttotlon to aubtraot 
tba oontenta of tba ragiatar wblob oontalna tba word > 1/2* Tbua, If 

* 1/2 ware atored In tba ragiatar wboaa addreae la 249 (or alaply ragiatar 

# 24^9 tba prograanad Inatnietion would ba au 249 , wbare jqi la tbe oodad 
fom of isytSeiisl^o 

Speclfleally tba Inatnietion lSLi4SL»aAn«» aubtraet the eententa 
of ragiatar ii 24^* Tbua Inatruetlona aay ba Ukanad to algabrale notatloar 
wbare tba eententa of a ragiatar najr ba ebaagad Juat aa aay tba waluaa aa* 
algnad to algebraic ayabola* 

People bave a natural but Inoorreet tendaney wban first begin- 
ning to prograa to write an inatmctlon aa an order plua tbe word to be 
operated on ratber tban aa an order plua tbe addraaa of tbla word} a-g*, 
tbe oorrect foni of tbe Inatruetlon above la au 2^^ , not an ♦ 1/2 * 
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Progruia ani orlglnallar written In atandard fom, that la, with 
tba alpbabatioal abbraviatlone of tba ordar ooda and ordinary daclMl ma- 
bara. aa in tba pracadlng axaippla* 

fiafore wmi can respond to the progran* the atandard fora mitt 
toiiahow be converted to binary fom* The taak of oonvereion it aesignad 
to the oompttter, inatmeted by a eodrvereion prograo vbicb baa been written 
onoe and for all* 

Tba prdgraaaer noraally doea all programBing In atandard fom, 
and need not concern biaeelf with the details of oonvertion« However^ be 
«4ll acre fully understand progrsaiiing if be is faailiar with tba font in 
wbiob nunbere end instructions are represented In VHJ» 

A 7.WX storage register consists of 16 binary digit positions t 
each digit posltioii nay contain a 1 or a 0} sea the drawing at the end of 
this sectione When a word within a register represents an ^na^'^'^tionp 
the first ^ digit positions are taken up by the binaiy ooded represents^ 
tion of the operation and the last 11 by the address of the word to be 
operated on* 

Vhen a word r<:f-rftff<;nt« a aygil^ftXt the first digit position, cal* 
led the sign digit. Indicates whether the number is poaitiwe (^) or nega- 
tive (-)• The laat 15 indicate its nagnitude* For a > nunber the sign 
digit is 0* For a * number the sign digit is 1. and the - nunber *s mag- 
nitude is represented by changing all 0*s and l*s of the positive nagni* 
tude of that nunber to I's and 0*s respectively* 
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Ttatts, * ? . it rcprtaenUd In • r»gi8t«r aet 
S 



•ign aagnltttde 

and the nagatlv® noabar -.X. *** 

8 



i|i|o|cii|iWi|il~xfflMiiLiJJl 



This represantatloB ie aora fully ascplalnad In Appandix Ao 
The nagatlva fora of tha nuabar is oallad tha ooaplaaant of 
the posltiva fom« and vioa varea* 

thus at any one tlaa a register oont&ine a vord ahleh ney be 
interpreted both ae a binary number and ae a coded inatructioa* \bhanever 
a word ie eent to the control element, it becoireF %n instruct ion« i^.ile 
vrheneTer it ie sent to the arithmetic element It ie treated ae a number* 
This will become clearer when exaaplef are considered later on* The fact 
th&t any rej^itter me^ cont&in either en instruction or a number interchange- 
ably makec for a more flexible machine* 

Km ie ao designed that arithmetic operatioaa are normally handled 
as if tha binaxy point of a musbar (oorrasponding to tha decimal point in 
tha HAoiaal aumbar system) is oonsidered to be fixed at the left of the 15 
digits which represent its magnitude* Thus only numbers of magnitude less 
than 1 are carried in a register* In other words, a register carrie;^ num- 
bers in the range between "*1 and ^1* 
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Of oooree VKI nuet be able to handle problena involving ma* 
bert greater than !• A nunber equal to or greater than 1 in nagnltude ie 
prcfultiplied bgr eone other numbei^ nomally called a eoale factor^ (nanally 
a poier of 2) eo ehoeen that the retult ie leas than !• Both the ecale* 
faetored nonber and the eoale factor itself are then needed to represent 
the original mabero The scale factor nay be assigned pemanently and re- 
nenbered bgr the programer, or it nay be stored In a separate register ^e 
the exponent of the power of ^ The progran nay be written to deal properly 
with the seals factor whenever it deals with the soale«>faotored nonber* 

Effective eoale-factoring of nunbers in the original progran 
and of the nonbers arising daring its execution is one of the progranBerOs 
more tedious tasks » If the nagnitude of a number resulting fron one of 
the arithnetic operations equals or exceeds 1 (called overflow ), the oonputer 
autonatically stops, gives an alara, and Indicates the point in the program 
where the overflow occurred* 

The binary nunbers representing addresses in Instruetlone are 
considered to be positive Integers with the binaxy point at the right-hand 
end of the word* 

When nunbers are represented in digital form, their aceura<9' 
is United by the nunber of digits used in the representationo The aoeuraoy 
of 15 binary digits (one part in 2^^ or 32, 768) is equivalent to that of 
about ^o 5 decimal digitso Greater aceurboy than this nay be obtained by 
using more than one register to contain a numbero 

An explanation of binary nunbers sufficient for general pzogran- 
ning purposes ie included in the Appendixo 
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Tb« four tetio •ItMott of tb« Uilrlviiid ooapator arts 
1* liipnt«ontr;ttt . to intfodueo prograM and data into tho 

eoffptttcr proper and to oxtraot roeultt of oalottlationt* 

2* SiSXMMMc 'or a program and data* (Stonge la now prlnarily 

oompotad of 304 regiatara; thia nnnbor la bolng Inoramaad 

to a fan tbonaond in tha naar Aitura*) 

3* Ar^^hBftio alanant f which addt, aubtraota* miltlplloay di- 

▼idaa whan ao Inatrueted* 
^* Control ,, which ooordinateK the ovfer^all perfornanee of WWI* 
Theee elenanta are inter-conneoted b^ a nain boa or eomunice' 
tlon link* The bus providee for the Interohange of Infovvation neeeeeary 
for liiWI*e autonatio operation* 

l^henever Ml is operating, the control elenent is generating 
a contlnuona and rapid anoceeaicn of eleotrical pulaea* Theae pulaea are 
fed to the iiqsut of an electronic dietribution ayaten through which they 
are aent to the circuita required to carry out eadh operation aa called 
for by the progran* 

The next four drawinga are illuatrationa of the baeic coaputer 
elesente* 
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Eo ttHAT yJWI DOES BITH PBDGRAliS 

KWI n^riMlly baa an input progran atorad in ita f irat 32 re- 
giatara* Thia input pro^rmm raada the progran to ba parformad into any 
otbar daairad ragiatare of tha atorega unit* At tba praaant tiaa, prograaa 
ara read into WI froa punohad papar tapa* In addition^ nagnatio tapa input 
end photographic f ila input ara baing davaloped* 

Aftar tha progran has baan raad into and atorad in VKI atoraga, 
it ia raady to ba parfomad* Upon a atart aignalg tba varioua inatruotiona 
in tha progran ara parfomad in tha aaquanoa in which tha progran liata than* 

Extamal data naadad during parfomanoa of a progran can ba 
fad to atoraga fren varioua typaa of aquipnant auoh aa radar aata^ tala- 
typa linaa, praaaura guagae^ atce 

Tha raaulta of a progran can ba raad out of ^I in irarioua 
foma* At praaant thay can ba typad, punched on papar tapa 9 and diaplayad 
on oaoilloacopaa* 

Tha Ordar Coda in i^ppandix £ givaa tha fUnotiona of tha ifarioua 
ordara* Tha ordara in tha fire t group aand controlling pulaaa to inp6t« 
output aquipnant* Thota in tha aaoond group provlda for tba transfer of 
infomatioa within tha ccarputar* Tha two ordara in tha third group nodify 
tha aaquantial opanttlon of oontrol to pamlt transfer of operatiolia fron 
one point in a progran to angr other point* All the orders in the last two 
groupa ara known as tha arithnatic orders* 

The following abbreviatlcns for three apeoial-purpose arithnetio 
registers ara used in the explanations given with the order ooda« £ach of theae 
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three registers hae a 16 digit vord capacity , Just ae do thoee in etorageo 
AC « aocunulator* The adding and shifting register in which 

SUBF 4^ products are a eonmulated and retained* 
BR « B«regi8ter* The shift 5jig register extending to the right 

of the aoouBulatui, used in the fonretion of products and 

quotients* 
AR « A«register« Primarily a buffer storage for words to be 

added or subtracted into AC* There is only one order, 

1§9 with which it is necessary to consider the contents 

of the A«>register« 

T«TJ was originally designed as an o^riaental Bodelo As it 
has been developicgp nore pemanest plans hawe been nade for it* NewerthO'^ 
less, it is still In an experimental 8tage« This is particularly true of the 
inpttt*output equipment and techniques for its useo For this reason the vem-^' 
porary orders listed in Appendix £, are being used at present « Instead of 
the first block of input-output orders listed in the Short Guide to Codingo 
The programs for particular problems will not be changed substaatiallyo if 
at all, by these temporary conditions* 
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Fo ^sssmim 

£v«Z7 tine a particular Job is prognunaed for ^111 tba pvogran 
boooBoa available to future progranmera, thereby eliainating iruch repe* 
titlon of effort* 

For inatance^ routine progr^na oalled aubroutinea, have been 
written for oonputing polynoniala, for coBputing aoat of the oomoa trig* 
onometrio Ainotiona, for oonputing the tquare root of a nunber, for arrant 
ging a eat of nuiibera in aaoending order of Bagnitude, for interpolating 
in a table of valuea, and for other taekso 

If a programer wanta to have one of theae taaka perfomed aa 
a part of a longer program, be need nerely oopy the aubroutine prograa al* 
rea^ available for the deaired taak, sake appropriate ehangea of addreea, 
and inolttde it in the longer prograa* 

An even nora ef f ioient way of handling aubroutinea ia bwing 
developed for IWI* It ia planned to oonpile a library of aubroutinea to 
be itored in a fom aooessible i>o ^fl* Any program requiring the perforaance 
of one of theae aubroutinea will inolude one Inatruotion directing it to 
the deaired aubroutine* The aubroutine will then be perfomed with the 
valuea aupplied to it by the aaln program* Upon oompletion of the aubrou^ 
tine 9 control will autoaatically be returned to the point at which it left 
the aain prograa* 

Eaoh aubroutine available to WVI in thia fom can be oonaidered' 
an extenaion of WKI^a baaic order oode^ ainoe the perfomanoe of the aub* 
routine requlrec btt op# inatruotion in the aain prograa* 
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The etraoture of « pzogras uy be clearly oharted on a flow 
diagrmm* A floe diagran eonslete of a aeries of atateHenta about what the 
progran doea* The atateaeiita are eneloaed In boaeea vhlch are oonnect^ by 
arroira whieh indloate the oourse of the progran* The atatenenta included 
in the boxea of a flow diagran nay repreaent one, a few, or naay inatrac- 
tiona, depending on the puxpoae at hand* 

Floa diagraaa are uaed prinarily for clarification jyurpoaeao 
Vihereae prograna are written in the language of the conputer^ flow diagrans 
need not be* They are helpf^ to the progranner in analyiing a problen 
and planning a progran* They are helpful to anyone intereated in what a 
progran doea* 
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PART II 
OEVELOPIfENT OF PROGRAMS FOR HWI 



Th« following oanplee have beon selected to illuetrate ho« 
prograat are developed for WWI and to euggeat other appllcationa and 
technique t* 

Fdr oonvenienoe of llluetration, Inatruotiona far prograne will 
be numbered as though they were to be stored eoheeoutively in regiaterag 
starting with the address lo Stored nunbers will begin at the address 
200 (an arbitrary designation)* Nunbers and Instniotions oan be assigned 
interchangeably to any address <> However, instructions DUSt always follow 
each other consecutively in the order in which they are to be perfomedp 
except when a special instruction within the prograa itself orders a de- 
viation froD this rule* 

When an order appears for the first tine in the following ex« 
anples the reader should refer to the Order Code in Appendix E for its ex* 
act descrption* 
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Ao ADDITION OF A SEUES OF NOVBERS 

Tb« first group of examples will illustrate the progrenning of 
arithmetic procedures* 

A»l-1* 1»ITHDUT SCALE FACTORING 

A program to add three numbere together is given hereo This 
program adde the numbers a, b, and c, and stores the result in a register* 
It then sende control to the first instruction of the next Job* It assumes 
that each of the numbers is less than 1 in magnitude, and also that at no 
time in the process will a number arise of value greater than 1 (i*a*y that 
no overflow will occur)* 

;»s^yVCl^^0B Mia^ 

1* ca 200 Clears AC; leaves a in AC 

2. ad 201 " (a+b) " " 

3» ad 202 »» (a+W-c)" " 

A* ts 203 " (a+bfc>" •• and In Register #203 

5* sp (address of register holding Send control to proper address 

* the first instruction of 

* the next Job to be done) 

Da^ 
200* a 
201* b 
202* c 

203* — - Contains (a-i'bt-c) after instruct 

tion in Register H is executed* 
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A.2o Mmi SCALE FACTORING 

llJuelXy« the exact values of a» b, and c to be used iritb a progran 
such as this are not known^ but their probable renge usually is known for 
noet practical applications* Such a situation could be handled by scale 
factoring to acconmodate the maximum possible veluee that could arise during 
the program o Scale factoring is done by multiplying a number by 1/2 enough 
times that the resultent (8cale-*factored) value is less than 1 in magnitude* 

For example, if it were known that a, b, and c each lie wihhln 
the range »5 to *5f since 5 x 2"^ « 5/8, they could be stored in the form 
of a X 2'*^« But the maximum possible value of the sum of a -** b -^ c Is 15* 
The largest scale-factored form of 15 which is lese than 1 is 15 x 2"^* 
Therefore, the appropsiate adaptation of the program to handle values in 
this range would bei 

1<=>5« (Instructions would be unaltered*) 

« 
t 

2CC. a X 2"^ 

201. b X 2*^ 

202. c X 2*^ 

203 i -»- (Contains (a+b+c) x 2"^ after instruction 
in register i'A is executed*) 

Corf^eepondingly, if the maximum value that could arise in this 
program were 20CC, an appropriate scale factor would be 2"^^^ because 
2000 X 2"*^^ « 1^ is lees than 1* 

IShenever scale factoring is z«sorted to, the scale factors must 
be remembered either by the programmer or in the program itself, so that 
final results can be reconverted to their actual values* 

A simple program like the above may be encountered as » «• 
part of a larger program, where different sections of the whole, o: 
program calculate and store a, b, and c, before execution of the in, 
tions which perfozm the additions 
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Bo EVALU ATION OF A POLIKOMXAL 

B.lo WSTHDIJT SCALE FACTORING 

The evaluation of the polynondal ax + bx + c for a particular 
value of X may be prograamed as follows: 

Instruction Effe^ 

lo ca 202 Cl^rs AC, a In AC 

2<> nr 201 ax in AC 



3. 


ad 


203 


ax + b in AG 


U* 


mr 


201 


ax^ + bxin AC 


5* 


ad 204 


iX^ + bx + c In AC 


6. 


ts 205 


2 
Transfer ax •»• bx + c to : 


7. 
t 

t 


sp 


(address of next job) 




1 








201. 


X 






202o 


a 






203. 


b 






20^. 


c 






205» 


oeo^ 




2 
Receives ax -^ bx ^ c 



This evaluation could have been accomplished l^ a different 8e<° 

2 
quence of operations! for example, by first computing and storing ax , then 

forming bx * c, and finally adding together these two values. However, 

this procedure requires several additional instructions and storage registers. 
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B*2. WITH SCALE FACTORING 

With the method Juet described for evaluating a pol^-s^mial^ an 
appropriately ecale^factored form of the first few terms of the Maclaurin 
series for sin Q , can be used to compute sin &• The first three terms are: 

where the angle ^ is expressed in radianso ^ is considered to lie in the 
range to -^ radians (O to 90®), but -|- is equal to 1»57« ^ is scale* 
factored to 9/2, whose value is always less than 1 in the range considered • 
Rewriting the polynomial in 6 as one in 0/2, there results: 

2(8/2) . -fiSZSlL ^ -22i£^SiL ;;j, ,1^ ^ 

For the limiting case of ^« •¥-, the value of sln^ is equal to 
1* Therefore 9 sin^ must be scale»faotored to ^ sinO^, and the final result- 
ing polynomial is: 

3 15 
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The 


corresponding 


Tnfitructiouf 


1« 


ca 2G0 


2* 


mr 200 


3. 


mr 202 


J,. 


eu 201 


5' 


mr 200 


6. 


mr 200 


70 


nr 200 


8. 


ad 200 


9c 


te 203 


10. 

1 

8 


ep next job 


1 

Ma 




200- 


e X 2"^ 


201. 


2/3 


202o 


2/15 


203. 


«..• 
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program to evaluate 1/2 ein & iaj 
(9/2 



- -|(©/2) + ^(<9/2)^ 

« -f (®/2)^ + ^iW^f 
- -|(©/2)^ + l|(e/2)5 

(^/2) - >|(®/2)^ + ii^/2)^ = 1/2 eina 



Receives 1/2 ein^ 

The oaxiinuB aceur&cy lost by neglecting the fourth and 
terms of the sin series is about oOO$ at 90°o For angles le«^8 
13° the accuracy lost is less than .00003 (or 2*^), less than ©an bs dis- 
cerned in a 15-diglt binary number. The inclusion of the fourth term ®f 

the series would decrease the evror of this ein P apj roximation to lese 

•15 • o 

than 2 for all angles up to 75 <> 

Any desired degree of additional accuracy could be obtained 

by using more terms of the series (making the equation inherently sbre 

accurate) I and more than one register to contain each number (making it 
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t)CS8'ible to fepreeent this greater accuracy In MI)o 

The raiiige of Taluee computed by tMe program could be e^teMed 
beyond 90^ by the addition of a fev^ orders to sense what q^iadrant the 
angle was in and accordingly to adjust the sign of the result* 



There are several other ways of finding polynomials to approxi-^ 
mate a function over a particular range (such e s the least squares method 
or Lagrange's method)* The particular method selected usually depends 
on the amount of storage required for the desired accuracy in the approxi-^^ 
matione 

Another way of finding values of functions is to have the prograia 
look up the desired value in a stored table of values, but this oonsa] 
requires considerably more storage space than does direct computation« 
The techniques for the use of stored tables is given in section He 



C" USE OF SUBRPpTINES 
C«l. TO FIND TAN &- 

If subroutines for finding 1/2 sin 9 and 1/2 cos O were already 
available to a main program, tan ^ could be found easily by using the re<^ 
lationship tan & « sin O/cob &■ o 

The program given here is valid for angles up to but tml including 
^50. for, without additional scale factoring, overflow would occuro A 
method for testing whether the angle & is actually less than A5^, end for 
using the program given here only if &■ is less than 15^ is given in seetione 

D«.l-2. 

Only the main program should be read at first* It indicates the 
over&ll procedure, with the sg orders directing the com|-;Uter to proceed 
to tx.e t^o subroutines and having the effect of leaving 1/2 cos ^ and 
1/2 sin & in the accumulator* The mechanics of the subroutine procedure 
is indicated following the program* 
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Eeport B-I56 

1. oa 200 

2. tp 150 

3* tt 201 

4. oa 200 

5. tp 100 
6* or 201 

7. •! 15 

8. tB 201 
tp next Jol) 



IfFBCT 



Main 

Prograai — — 



0/2 



So 1/2 GOB eu'broutln«, 0/2 Btill in AO. 

Tranifert 1/2 oob to 201. 

0/2 

To 1/2 sin subroutine, 0/2 still in AC. 
^/§g|gg - tan in B-recister, 

tan in AO* 

Transfers value of tan to 201. 



1/2 tin 
subroutine 



( 100. ta 110 

( 

( 101, ts HI 

( 102, mx 111 

( ' 

( • 

( ' 

•—..•4 110. sp _ _ _ 



Puts return address of nain prograa in 110. 
Transfers 0/2 (still in AC) to 111. 



1/2 Qos 

iubroutine 

( 
Itorage for ( 
■ain prograaiC 



( 

( 111. — 

( 

( 112. 2/3 

( 113. 2/13 I 

( lli^, — ) 
I 

I 

130. ta 160 ) 

• )-— instruotioas 

• I 
160 • sp -— ) 

I ) data 

• ) 



This Bueb of the subroutine ii identical 
«•-> with previously developed progran for 
1/2 sin 0. 



SaBe procedure as in 1/2 sinO 
subroutine. 



200. 0/2 

201. 



storage for values of 1/2 cos and tan 0. 
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The <g2 ^'^^ Ml orders are introduced in the main programo When» 
ever the ^ (divide) order ie need, the quotient of the operation is left 
in the B-regiater« ^ is usually followed by el 1^ (shift everything in 
the B-register and in the AC^15 digit positions to the left), which places 
the quotient in the AC* 

^^ 8p 150 In register ^2 interrupts the sequential perforiiBftiVJe 
of inatructions in the aain progren, sending control to 150* During the 
performance of the sp 150 in register #2, the next address 3^ following 
the jgs order in the main progrem, is placed in the A-register* The ta 160 
order in 150 (next to be performed) takes this address 3 from the A-»register 
and places it in the address section of the instruction 160, making it sp g> 
Thus, upon completion of the 1/2 cos & subroutine, control is returned to 
3 in the main program with the value of 1/2 cosS- in the AC* 

Similarly sp lOQ in 5 sends control to the 1/2 sin^ subroutine; 
and the combined action of Iba IK) in 100 and S£^ in 110 returns control 
to 6, with the value of 1/2 sin V in the AC 



When appropriate subroutines are available, as was supposed xa 
this example, the prbgrammer's job is greatly simplified* His attention 
is confined primarily to the main program in which he determines the proper 
use of the subroutines* 

A given subroutine may be used as maqy times as needed in a main 
program* Each time, an ££ instruction in the main program to the ta in- 
struction beginning the subroutine is all that ie required* 
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D- CONDITIONAL PROGRAMMING. FLOW DIAGRA&S 

The BubproEramning order JSB makes possible versatility In program- 
ming lay Interrupting tbe sequential performance of Instructlone* The 
conditional program order SU provides additional versatility by letting 
the computer decide whether or not to Interrupt the sequential performance 
of instructions* 

If the nufber In the AC when a sSi instruction is being performed 
is negative, the SESL operates exactly as an jb* If the number in the AC 
is positive, the 52B 1* ignored and the sequential performance of Instruct 
tions is continued* 

Thus, 'VTiI chooses between two courses In a program each tine a 
jSB inetructlon is given* It Is in this sense o and only in this sense, 
that VIWI makes decisions* 
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Peee 'i2 



D.l. USE OF CP ORgER F OR DETERBINING IF Q'lS < L^ 



As ft typieel exiiplft of tk« um of the jQ2, ooneldftr the proeedlng 
proerftv nhitfb comptttfte t&n &for valuos of O lost thftn irA rftdiana U5^)» 
Aeeuno that Taluos of O greater than 45^ are also eoBetiisee available in 
reglater 200 (in the eoale*'faotored fom ^/l^ before the operation of the 
progrcn* The following three inetruotione inserted ahead of the nain pro- 
gran would have the following ef feoto For ^ less than v A "^b** program 
would oontSnue with the evaluation of tun & • For values of ^ equal to or 
greater than v/A, control would be sent to tone other progran* 



ca !^ H^^ * ^/^ 



2« su 200 

3 



^ (address of next job) 



ir A « 1/2 

1/2 (irA-(9) 

Goes to next job onl^ if (ir/i ^Q ) 
is negative, otherwise oontinues with 
tan & prograa* 



'* The abbreviation RC should be read "Address of Register Contclsaing* 



D«2« FLOff DIAGRAM FOR TAH 9 PRD;31AM 



liith this addition, the tan & progran could be represented on a 
flow diagrasi ass 



Is e <^5®? 



VSL 



yes 
^OBpute and Store 

1/2 cos e 

[CoBJrute i/2 sinBl 

Compute and Store 

j-/ft i ^1^ , . ^ . ai tan& 
1/2 oos ©• 



.to next job 



1 



to next job 



Eaoh tine arrows narked "yes" and "no* leave a box of a flow die* 
gran, the usf of a £B in the program is indicated* 
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D.3* GENERAL REMARKS ON PSB OF CP ORDER 

Th« ssk order uaed alor^ ehooses between two altematiire oouriee 
00 tbe baele of the eign of the nunber In the aooamulator* The previoua 
esavple need the soi ^o Mke this oholoe on the baaia of the aign of the 
differenoe between two noaibera* Thua, the fl2,uaed with other ordara.aade 
the ohoioe on the baala of idiioh of two nunbera waa larger* . With other 
prooedurea ualng the SB i^ i* poaalble to irake ^arloua aeleetlona audh aa 
the largeat of a aet of nuabera* With the SEBt order oheoking the eontenta 
of a ooonter it ia poaalble to progran the cytilieal repetition of a aet 
of ordera a predetermined nunber of tinea* 

The following exaqplee aenonatrate auoh teohf.iqtte8* 

The sk order nay be used In a Tarl^tv- or other waja to diatingttiah 
between different aituationa or thinga which have nnnerioal valuea aaalgned 
to than* Sueh thinga aa problena of logic and varioue economio ganea nay 
be developed bgr having the sg^ diatingniah between the eigne teaeooiated with 
tme and falae atatenenta and different eonditiona* 
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UA* &mx^m Qf m imu^ py a ^p* qf ^m^ * d^bct ygrrop 



Thle progras finds tbe largttt of three annbere, •» ^t a^ o 
end etoret it in register 203* 



-A- J^ .^ 


1. 


oa 200 a 


a 




2. 


•u 201 (a-b) 


(a-b) 




3. 


qp 7 po». 
ca 200 a 


neg. 

1 
! 


pot. or 
neg* 


5- 


ts 203 a 






6. 


ap 9 




i 




?• 


ca 201 




b 




6. 


t« 203 

> 


1 


b 




9. 

10. 


eu 202 (a-c) 
cp 12 pos. 


(b-c) 
pop. 


,(a-c) 

or(b-c) 

neg. 


11. 


^ next job 
■•P - -. - - 






12. 


ca 202 




c 


13<. 


te 203 




c 


u. 

• 


.^ next job 
ep -. 






200. 


a 






201. 


b 






202. 









203 • 






Contains «, b, 


or c efte 



1-3 



Floe diagreft 



1-3 



^-6 



7-F 

9*-lC 

11 
12~13 



Is a > b? 




yes 


■ 
no 


i 




Store § in 
203 






1 










y 




Store ji 

in 203 




■ ■ •■■ ■■ •) 





Is OK* 203 > c? 


yes 


no 

f 



[ To next job 
]~Store 5 



ITTT^ 



i 

To next job 

^he abt^reviation OR 
should be read *'con* 
tents of register ^ ^ 
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This progran oould ba axtandad to find tha largatt of any nuobaF 
of Taluaa Iby tha addition of a aariae of ooqpariaon ordara for aaoh addi- 
tional yalua* Tha Baeaaaary rapatltion of idantiOfiX ordara (vith diffarant 
addrase aaetiona) in auoh a progren ia avoidad in programa aaing ooontara** 
to ba daaerlbad In aaotiona D«6«-7* 



Tha nuBoar ie rapraeantad in tiFT aa aithar ^ or aa « 0« 
Noraally a in a prograa raf ara to ■*' C* Bbiiavar« «han ia tha raault 
of an addition or a aubtraotion ordar^ the raault ia laft in tha aoounu- 
lator in tha fora * 0* Tha sk ordar ia aotuatad by • aa by any othar 
nagativa nunbar* thia point mat ba conaidarad ahanavar tha £b ordar ia 
uead in ehaeking ralative nagnltudaa of mobara, aa in tha praoading %t* 
ABpla, and in oountara, uaad in axaapiaa to follow* 

Tha axoaptiona to thia rula ara that <-*' 0) plua {* 0} laat — 
(♦ 0) in tha aeouaulator, and (^ 0) ainua (• C) alao laavaa (■»• 0) in 
tha aoeonulator* 

In tna praoading asaiipla, sif « and ^ vara aqiial in valua and 1 
largar than stt tha path followed in tha prograa would ba tha one indioatad 
for ^* Thia ia bacauae a <» b would equal •« 0, which would actuate the en 
Inat ruction in regiater H3* 
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D.6. OTyyiMi Q TO Pi W -fi tff C 



T.1^I can be programed to count up to a given number j^* Tbeee two 
prograne are called ootlntere* It ie neoeeeary to etore (n-1) in register 
jr201 to Bake the ojole be perfomed jq tinee* because of ^O* 



1» AC Oft K'th gyrty 



1« 


Cfc 200 


2* 


te 202 


3. 


ao 202 


/►• 


eu 201 


5* 


cp 3 


6. 
t 


ep next Job 


2C0. 


C 


201. 


(n-l)(x 2*^^) 


202* 


k(2"^5), counter 



k - (n-1) 



1-2 



L»« 



ST^ 



react counter to 



add one (1 x 2*^^) to kl 



r: 



Xe k ae large as n; 
I.e., ie k > (n-1)^ 



no 
to neacV 



yee 
job 



k -> number of tinee the cycle hae been 
repeated 



Eeeentielly thie propraii counts from to no The following program 
counte from «n to 0* 
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!• 08 200 
2* it 201 
3* to 201 

4« cp 3 

•p aMrt Job 



5« 
I 

t 
200< 



to AP 9ft )^'1^l> oyglt 

-(n-l) 



-(n-l) * k 3 

Mg* until k « n 



201. .^r: 



(ii-l)(2"^5) 

-(n^l) ♦ k (2*'^^), counter 



START 
~1 



1*2 reset counter to •*£ 



add one to ( k*n) 1 

Bae tbe counter I 
reaofaed 0? | 



no 



yea 



to next job 



The second of these two oountere requires two fever registers 
than the f irate Therefore it would be used in preference to the first eac» 
oept when it nay be desirable to keep track of the value of ]| specifically* 



These prograns can be adapted to repeat a certain process jg tines 
by having the counter cycle Include the desired proceaa« 
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C«»7<»' CYCLICAL PBPGRAM FOR FINDING LARGEST OF A SET OF MtmBERS i 

CHANGpG INSTRDCTIGNS DDRING OPERATION OF A CYCLICAL PF^GiUlfl 



Thit program uses a counter in finding the largest of a oonbere 
atored in IC^I and abova how nwi can change its own 1 net root ione* The 
counter providea for the repetition' of a aeriea of conpariaon orders until 
all J2 nunbers have been checked* The '^^neral fom of the progran nay be 
planned on a flow diagranc 

Start 

L 

{Reset counter* 



[Prepare to coaipare the lat two n unbere 



First tise thru cycle ^ 
compare the lat two 
numbers and store the 
larger« On successive 
triale^ compare the 
largest number yet 
encountered with the 
next unchecked number, 
an d ato re the largere 



Add 1 to co unter* | 



JDoes counter indicete 
that ell n numbers 
have been compared? 



no yes 

to ncart job 

This procedure suggeete itself because it ie similar to that fol- 
lowed by a person selecting the largest of a series of numbers o 

The tjj order le introduced bei-e* B«th the td order and th© ^ 
order are used to change the address section of an Instruction so that the 
same order may be used to operate. on the contents of different registers 
on different runs through the cycle o 

, * xP'® possibility of programmirs? mi to change its own instnjstloo? 

16, as is the: use of the 5J2, one of Ite rrc'st important f*-ator-#s« 
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n 

3 

5 
6 
7 
8 

9 

10 

U 

12 

13 

U 

*15 

*16 

*17 




cp e 

^ next Job. 



Cyolio«l PxognuB For 
Finding Lctfott of a Nuaboro 

Botet oountor 

»i 

Store Jt. at tbe largest nurber jet found* 
iddref B of X2 

Transfer addreee of x^ into tbe oa order In rege 8 

II m mm nniiii ***11 

Comparieon orders 

Store newly found naxlnua 

Prepare to ooopare next noiiber 

Increae« counter bgr 1 

Have all jq numbers been cheeked 



X sax* 



Data 

*200 

*2C1 
202 
203o ri 205 
204<> Xj^ 

2050 x^ 
t * 
I 
I 

203 * n« X. 



(n-2) X 2"^^ 
counter 



) 



Storage for largest nuiiber yet encountered 

Addreef of x^* Tbe binaxy coded fon of the order 
ri is OCOOO. 

n nuBberSy x. through x 



^Registers used for the counter in thle prograne 



Inetruetione 1« 2, 15» 16 ^ and 17 are the eane as the instructions 
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@f the -n to counter progran given earllero The counter reset instruetione 
BTuet be perfortted outtlde the oaln egrcle, while the eo, cp, and ep inatnietione 
BTUBt be perfonted within the nain cycle if they are to control the nonber of 
tinea the sain cycle it performed* 

Whenever a counter ie used it ia inportant to oonaider the exact 
nunber of tinea the controlled cycle ia to be perfomed* To oonpare d nuabere 
ae in this program, the conparieon mat be perfoxved (fi-1) tioiea* This requires 
that the counter here be reset to ^dt-S)* If the counter instructions within 
the nain cycle iiad preceded, inatead of followed, the eoopariaon instructions, 
the counter wo&d have been react to «'(ir^) f^^ Vfi*l) perfomanoea of the con* 
parison inatruotions* 
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D.S. FtOli DIAOLil FO B ARRAMGIKG A SET OF NUMBERS IN ORDER 0¥ MAGNmPE 

Tba prtceding program may be extended to rearrange a aet ox £ 
nombere in detcendlng order* One way of doing thle ie to find the largest 
number in the aet and erohange it with the firtt number in the set* Next find 
the largest of the set of nrl remaining numbers and exchange it with the 
first number in the new set* Continue the process until all numbers have 
been arranged* 

A program following suoh a procedure would oonsist of two main 
cyoles, one within another* The inner oyole would be the same at the main 
cycle of the previous program for finding the largest number in a set* The 
outer cycle would in addition arrange for the interchange of numbers and the 
elimination of numbers already arranged in order 6f magnitude from the next 
set to be checked* 



A flow diagram for this process may look as follows s 

START 



Prepar&tion for prograir for 
finding largest of a set of n um bers * 



r 



: 



Progr&m for finding largest 
of a gft oif BUffihirt* 



Have all nuabere in set been checked? 



no 



1 



yes 



Interchange the largest number with 
the first number in the set* 



le set by 1 swmber hf eliminating 
( the previously found largest no mber* 

JI.IUHIPI - - - ^ 

' Have all numbers been arranged in order of 
I magnitude, i«e»» have n»l sets been checked? 



no 



J 



yes 
to next Job 
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The Interchanging proceee requires ^ and ^ and j^ orders* 
The jgfi takes the inetruotlon which contains the address of the register of 
the largest nunber, the tg transfers this address to a ^ _ . . instruction* 
This jifi follows a cft RC (first nninber in the set )* Thus the first nuuber in 
the set is transferred to the register that originally contained the largest 
number in the 8et« Sinilarly the largest number is transferred to the register 
which originally oontained the first nunber* 

The following sets of instructions added to the previous pro- 
gran would provide f6r this Interchanges 

First number to original Largest nunber to original 

register of largest^ nanbgi: pggMtf r gf f iret n^jff^r 

n« ca EC address of largest n'* ca 3 

number 
n ^- 1* td n + 3 n* ^ 1* td n* + 3 

n + 2* ct 20^ fj' + 2o ca 202 

n ^ 3» ts ^ . ^ „ ft' ■*" 3» ts\^ _. _ ^ 

To eliainaiethe first nunber from euoceesive sets after the in** 
terchange, it le necessEry to give &n ^g order to :>.n8tractlon 3> and to 203 1 vhich 
set up the first register of the new set to be cba ikedc 

The addition of one extra counter will check to see if all n » 1 
sets have been checked* 
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It l8 frtquently ntMtiary or oonvoniont to •olvo nonorloal 
pvobloBf bor itomtlTo mtttaods— that la, by mtboda whleb vako ittooaativo 
trial attanpta at a lolutlon until a aatiif aotory aolutloa it raaohad* 
1iif:i*a rapid operation aakae iuob aatboda practioal* 

&.i. TO \mS%n LiWAM yymtfiro wq<^es§to triam 

A iiaplt iterative prograin for finding the equare root of a poai- 
tive xnuBber x leea than 1 ia given here* Beoauae oa 'jiirltttt a thia.pro*> 
oedure requiret aeveral thoueand trlala to find the acuare ropt whieh oan 
be evaluated readily in other wagra, it ie of little praetieal iaportanee, 
but aervea at an lUuatration of what oan be^done bgr iterative Mttaoda* 
Let '/F Indicate euoceaaive trial valuee of Vx* 



^l&a^p^etiona 




lo 


oa 202 






te 2C0 




3- 


ao 200 




il- 


nr 200 




5- 


•u 2C1 




6. 


cp 3 




7. 
t 

1 


ap next 


Job 


1 







x« 

x* - X negative until x* becoinee 
larger than x 



200. L'^ Initially 0, finally \'% to within 2^'^^ 
201« X 

202* 

^be ao 2pQ inoreaeee 1^ by tL^ on each eucceeeive cycle until 

x« juat exceeds x« The flret value of fi^. tried ie 2"^, and x» here le 

'•3C 
'Si ^ which would appear as > in th« accutmlator* In fact, x' would 

appear ae ■^ in the accumulator untfi j'?^ « 2 ^^ wae tried, giving 
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x« « 2 "^^ Thus, 11 the above progran were ueedp 2 "^^ instead of 0, 
would, be stored in register f^202 in order to conserve computer tine* 

xnie progran nay be used for values of x greater than 1 by scale 
factoring x by 2"^* This would leave x(2)*'*/* in register #200« 

If 2~^ were, snaller than the accuracy required for a particular 
problen, the flfiJgfiQ could be replaced by the three orders o^ 20^ ,. ad RC- 
Cdeeired accuraov) . and %p gpQ. Alternativexy, an sx n and a ts 200 

following the ao 200 would give a* difference of 2 x 2" between successive 
values of Yx^. 



With iterative programs it is important to consider the amount 
of computer operating time required to execute the program* An average 

operating time of IC seconds should be allowed for the execution of each 
instruction* This is equivalent to the performance of 10, COO instiuctions 
per second* 

The irSTprogram given here would require up to about 12 seconds 

for values of x near 1* This estimate is made on the basis of about 2^ 
or about 32, COO iterations of the A instructions ^ 3-6* If between suo» 

cessive trials steps larger than 2 "^ were used^ this time would be out down 
proportionally, at the expense of accuracy* 

The next program for finding the /iTby Newton "s iterative method 
converges on the solution very rapidly* On the other hand it requires 
about twice as much storage space* In many problems where alternative 
methods of solution are available the choice will depend largely on opera- 
ting time versus storage space* 



These programs of iterative processes are given. as examples of the 

general procedure used with such processes* T* e addition of (n-2) Instruct 

tiona of the form mr RC /p to the program just given would give a program 
to find the nth root of x* »Slirilarly, Newton «s method, used to find a square 
root in the following program, may be developed to find a root of any nth 
degree f(x)o 
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E^S. fx BY NEWTON'S liffETHOD 

Newton* 8 nethod for finding e root of an equation f (x) ~ utee 
the fomala; 

where x^ ii the nth approximation to the roo*^^ and f*(x ) it the elope of 

the curve at x • Succeeaive valuea of x. are found at the interaection of 
n n 

the alope f <> (x^) with the x axis* 




^•<vi) 



The aketeh indicatea that succeseive approximationa converge 
rapidly on the value of the rooto Therefore a high degree of ecouraey can 
be obtained with a few iterations. 

To find the equare root ,^ a nunber % let f{y.)^x •a«0« The 
formula for aucceeaive approximationa to the positive root of -this equation 
becomes s 

This program aasumea that 2"^ < a < 1 - 2"^« Starting with an 

initial X- « 1 - 2°^^^ it finds successive values of x^ until |x^ - x ., I <2"^« 
A n I n ff^i I 

The last value found for x^ will be an approximation to »T, aocurete to within 

Thua the program must have stored initially three valuea— the 

-15 
number n, the first approximation x. » 1 « 2 , and the. desired aceuracy 

2 "^o There ahould be reserved a register to contain s * Another register 

is necessary to contain the partial result, 1/2 (* " 3f-^ 

*n '• 
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Tji^p^ftiotion 


le 


oa 201 


2. 


ts 203 


^* 


ce 200 


k* 


dv 203 


5. 


8l 15 


6. 


sn 203 


7. 


sr 1 


8. 


t8 204 


9- 


ad 203 


10. 


te 203 


n« 


OB 204 


12« 


an 202 


13* 


op nesEt job 


1 



Q 


8P3 


200* 


a 


2OI0 


Xj « 1 - 2*'-^ 


202o 


2-U 


203« 


— •» Ml 


204. 








a/x ' in B.-regieter 

a A In AC 
n 

l/2(a/x^ - x«) « X 



n*l 



Vl 



'nn 



Return for oaloulation of next x_ 



Uaed for x^^, finally [a 

Used for l/2(* - x^) 
n 
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Fo GENERAL NATURE OF PROGRAMS lOR SOLVING £IMULATANEX)US 
EQUATION AND FOR INTEGRATION AND DIFFERENTyATXON 

iol. SIMDLATANEOUS ALGEBRAIC EQUATIONS 

Simultaneoue algebraic equations nay be solved on WWI 1^ a variety 
of metbodfio The investigation of such nethods is a problem in numerieal 
analysis^ and no attempt wil3 be made here %6 consider them since they 
represent essentially a mathematical analysis problem rather than a program^ 
ming problem* Onee an appropriate method is selected for a problem ^ the 
programming consists of scale factoring to avoid overflow and expressing 
the process by a sequence of instructions « However, the choice of method 
must take into considercrtion such thiigs as whether the method will always 
converge on the solution^ the aceuoulation of round'^off error, required 
storage space, and operating time* 

Simply to suggest the ilature of the problem, consider the se^ 
of simultaneous linear equations s 

Several methods of S'^lution suggest themselves • A direct sub^ 
stitution method would be to solve the first equation for y « f (x, a,, b„ , c«)o 
Substitute this esepression for y in the second equation and solve it for 
X * ^^^l* ^i» *^i» ®2» ^2» **2^' ^^^^ * similar expression for y- Then 
the program would oonsist simply of the direct evaluation of these express 
eions for x and y in terms of a, b, c« This method resolves itself into 
the same process as the evaluation of x and y l>y the irethod of determinants 
or even by more general matrix manipulation methods* In some problems these 
latter methods may be appropriate for WHQ* 

A less direct method of solution is the iterative one using the 
relations ; 

*1 

where the si^bscript i refers to the ith successive trial for the solution 
of X end yo Initial trial values x. and y. are assigned to x. and y., 

from which are found x^, y^, x^*** x „ y successivelyo The iterative 
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prooedure is itopped when both of the differeneee x^^^^ "^ ^1 ^°^ 

I ^1-*-! ^ ^i I ^^°'® ^^'^ than the Daximin allowable error^ 0° This nay 

be done by applying the sSl to \x^^^ - x^| • e, and to jy^^^ * y^j - e* 

START 



^ 



Transfer initial trial ▼aluea of 
X. 9 ^1 ^ regiatere for x., y.* 



Calciulalo X, 



i+1 



\ ' Vi 



Find and store X4«.* - x.« 



Store x.^. for the next cyele* 



I 



Celottlate 



'in 



*'2'*'g|4.i 



Find and store y-^. - y.« 



I 



&tore y^^. for the next oyole« 



pnrTx 



;a3 



m. 



VL 



"eT 



mr 



no 



ym^yjO •' 



yesj Tno 
sp to mat job 



Iterative aethods for solving simultaneous equations requin 
longer oonputer operating tinOp but nay be easier to set up and progxwi» 
and in sone oases provide the only praetioal neans for solving large wy^- 
teas without undue round«of f error* 
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F.2. imm^^rm 

Integration and differentiation nay both be approxinated in VWI 
by replacing these oontinuoue proceeees by discrete diecontinuoua 
0U8 numerical procedures—the only kind ^t purely digital deyices can 
handle* 

Integration nay be aprroxinated in ^VI by various processes of 
sunaation* For instance, the area of the figure below oould be approzinated 
by the standard methods of summing narrow rectangular elements of the area, 
by the trapeaoidal rule, or by Simpson's parabolic rule* These methods 
for approximating integration inwolwe a series of steps, each of which may 
be handled by tivi; i*e«, the evaluation of f(z), multiplication to find 
the area of a finite strip, and addition for the summation of these ele- 
ments* 




-> X 



^jt, 



The diagram indicates the approximation made by summing narrow 

D 

rectangular strips* Here the area between C and D is^ f (x) ^ x* The 

C 
f (x) dx is represented 



difference between this and the actual area 



I 



by the area of the small segments between the top of the rectangles and the 
curve of f (x)* This error, called truncation erroy . is inherent in all 
discrete methods for apparoxivating continuous processes* The truncation 
error may be reduced by reducing the discrete interval A x or by using 
more accurate approximation formulae* The former usually increases com- 
puting time and round-off errors; the latter requires more storage and is 
harder to program. 

Cyclical programs seem most natural for evaluating integrals 
on mil* The flow diagram shows the general prooedure for a program to 
find the area under a curve f (x) between the limits x « C to £• Let the 
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total accunulated area be represented by^AA., where AA. le the incremental 
elements of area taken at the value of x equal to x.* 



Start 
J, 






Set Initial value of 
X. to its lower limit 


of 


C 


w 


Set^AA^ to 







Evaluate f(x.) 



Evaluate A A. by an 

appropriate formula 
aeing x. and f(x.)* 



I 



Add AA^ to ^AA^« 



increase x. ty A x 
for next cycle* 



Te x^ > D? 



I 



no 



yes 



sp to next Job 



2aa, 



The result ^ A A. will be left in the register used to contain 



Uore complicated integrations wil]^ call for programs of this 
general form* 
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FOo DIFFERENTIATION 



rifferential equatione, like integration, describe oontinuoue 
proceeees which must be traneformed to forme which can be treated by discx^te 
prooeesee* Among these are difference equations, extrapolation formulae, 
and power aeries expansionae The selection of appropriate methods is more 
a matter of mathematical analysis than of programming, and again only a 
relatively simple but important exaaiple will be given in the next section* 

Simultaneous differential equations are handled by first trans<» 
forming them to simultaneous algebraic equations* 

A program to solve a set of simultaneous differential equations 
by the use of linear extrapolation formulae is given in the next sectiono 
In general, extrapolation formulae carry an approximation to a curve along 
short segments of lines parallel to the curve which Is being approximated* 



Consider that the first derivative 



(-g-) 



at a point on the 

curve of a given function f (x)— defined by differential equations- 
is the value of the slope of the tangent to the curve at that point* The 
slope at this point may also be expressed as the ratio of finite lucre* 

ments (^;^)£* These finite increments are used to replace the different 

■tlals in the differential equations* 



y + 






/i 



ext •>"/U^ef C^-'*'*. 







X 



y^^2 * ^i * ^dx" ^1^*^ ^~^^l " ®^op® f o«n<J ty evaluating 

derivative of f (x) 



'W 



x^ "»■ Ax 



; A X -s constant 
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Go SIMOUTION AND DISPLAY 

Gel. SOLUTION OF DIFFERENTIAL EQUATIONS DESCRIBING MOTION 
AP O^Cim^CO PE J3}SPm OF THE fATH 

Thle example of the handling of a set of differential equations 
will aleo be used to illustrate simulation of a physical proeest hj vn* 
Use of the output display orders for controlling an oscilloeoope will be 
introduced* 

The differential equations of a projectile, or any moving body 
in space, acted on by the force of gravity only are: 

-fi-1- « - g (acceleration of gravity) 

lit~ ^ ^ (constant horizontal velocity component 

Linear extrapolation fonmlae used to find successive points 
(x^, y^) of the path of motion at tines t. are; 



Xo^^ -X. +v At*x position component 

^1+1 ~ ^1 ^ ^\^i+l^^ '^ ^' P®''^*^^" cotrpenent 

(v )^^j^ - (\j^ - g^t » vert leal velocity cotfiponent 

The program solves these equations for 3c-^», y.^„ an«^ displays 

the corresponding points on an oeeilloecope» It further displays a hori- 
zontal axis at y « and arranges to have the path of motion r«*present that 
of a body bouncing on the horizontal axis» The coefficient of r^stlUition 

H is equaJ to - \f^ Wo^*^^^1 The a 1^ "«d to d,t,«,ln,e .hen 



the value ©f /, becomeif^ negativeo Ih^n the fwrrerpondiag doira*ard value 

of (v, ), <• we?t^rli«c by h tc7 g«t. (v } which will be direet«t- iu;;««rc3 
since R if negative* a sketch of the display follows th« prograaio 
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BpTTNCTWn B/ILL DISPLAY 



1. 


ca 


207 


2* 


t8 


201 


3» 


oa 


208 


A- 


te 202 


5. 


ea 


209 


6. 


tfi 203 


7. 


oa 


201 


8. 


ad 210 


9. 


qh 


201 


10. 


ca 


*8;12 


!!• 


qd 


212 


12<. 


ee 


202 


13. 


cp 


17 


u. 


c« 


203 


15» 


mr 


204 


16« 


t» 203 


17. 


ca 


203 


18. 


8U 


205 


19. 


ts 2C3 


20. 


mr 


206 


21. 


ad 


202 


22. 


qd 


202 


23. 


eifi 


201 


U. 


8U 


211 


55. 


cp 


5 


26. 
1 

1 


•P 


1 


201. 


.*!. 


202. 


^i 


203 « 


•7 •«-"«• 


20i^. 


R 




205. 


g^t 


206. 


At 


207. 


^1 


I 


20«.. 


^I 


s 

f 


209. 


(v 


). { 




y ij 


210o 


\' 


At 


211. 


^U] 


J 


212. 


+ 






'1+1 



*1 ' 'x^* 



Set horisontal deflection of oscilloecope to x^^n > 
transfer x^^^ to register #201. 

+ 

Set y deflection to ■♦• 0, display point (Xj^.!, ^ O) to 
form a point on the horisontal axis, 
y . ^ Determine if last point calculated lies above or 
/ below X axis and proceed accordingly 

<V^b 
R(v ), 

y 

y 1 



,. (v ) « new (v ). directed upward 
b y'B. y'i *^ 



g/it 



Display point of path at (x-^-, ^±^1' 

N 

^^ I* ll'<.^L» '^^'^^'^ ^ calculate next point on path 
Return to repeat diso3«y 



) 






Initial condltione 



yppt-r Llirtlt of X. 1 greater tha;?ii x^ In magnityd«; 
y value of horizontal axis 
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Of course aetual nunerieal valuee for the stored constants must 
be used in the final program given to f:wi« These values would be scale- 
factored and proportional to those of the siimlated process* 



A program such as this could include fomulas for modifying the 
initial oonditions of position and velocity so as to control the actual 
path of notion during operation of the program* For instance a relationship 
between the final conditions at x^ and the initial conditions could be pro» 

grammed, which would serve to adjust the initial oonditions to meet certain 
desired final oonditions at x^. « Further, it would be possible to have the 

program provide for the variation of x., ^^^h externally introduced datae 

In this way the program could be considered to both simulate and control 
a physical process* 

At any rate, the parabolas displayed as successive points on the 
oscilloscope screen represent solutions of the original differential equations 
for different sets of initial conditions* 







Results of programs may be displayed In a variety of other ways- 
At present results may be typed on a typewriter or punched on paper tepe* 
These latter two forms of output require special conversion subroutines 
which are available* 
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H. USE OF STORED TABLES 
R.l. GgMPM> PPQPEPypg 

The general procedure for looklag up a value etored in tabular 

fom In V<VX storage is given here* Consider that a stored table of values 

consists of A values which are functions of x, x ranging f roa 1 to n* The 

progran begins at a point where a value of x equal to x. is already in the 

accumulator* The last o^tier leaves the corresponding functloug f(x.), in 

•15 
the aoouaiulator* The values of x. are scale factored to x. x 2 ^^^ 



Ins-^IVctlOfte 
t 


£f£j?SJf 


Ms 

t 
1 


I 


t 
U* 




Xj in AC 


215. 


ri 215 


15* 


ad 215 


ai + 215 


216. 


f(Xj«l) 


16* 


td 17 




217. 


f(x^='2) 


17. 
« 

t 

t 


ca ^ 


f (x^) in AC 


218. 
1 
< 

215+ n* 


fCx^-n) 



Notice in addition to the table of values that it is necessary 
to store also an address to indicate the beginning of the table » 

Normally it requires less storage space to evaluate an analytic 
function than it does to store a table of its v&luee unless only a few 
values are needed* On the other hand, the use of tables nay save consider^^ 
able operating tine, which may be needed in certain control applications. 
A method for interpolating between values In a stored table of analytic 
functions is given in the next section. 

Tables must be used for looking up non^analytic functions of x, 
i*e*, values or words which are identified by the value of x but which cannot 
be evaluated from x in an equation* For example, such tables are used for 
conversion programs to get from one coding system to another* It is with 
the use of suc^ a table that l^I converts decimal and binary characters to 
binary when a program is read into storage* The same conversion tables 
are used to convert results back to decimal and alphabetical characters 
when results are being reed out of the computer* Such tables are also 
needed in al^ applications involving accounting, statistical analysis, 
correlation of data, etce 
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H»2. LINEAR INTERPOLATION 



The fonmila for finding an unknown function of x, f (x), which 
liee between two known functions of x, f(x^) and fCxj+j). by linear in- 
terpolation its 

f (x) « f(x^) + D ^^^in^ - f<»i)J 



X - X, 



where 



w =* 



*i+l " '^l 



X - X, 



Ax, 



ie 



the fraction of the interval Z^x^ that is covered by the Interval (x-x^) 



I*- 







j / ' i * 
I iJLl 



^r^ K*> 



■\ 



VrrvAA;, 



X 



\. 



_i • ( ♦- 



Tu the program given here using this interpolation formula, values 
of X, are considered to be positive integers* The value of x « x. *■ mAx. 
for which f (x) is to be obtained hy interoolatlon is represented 1n p reg"- 
Ister as X X 2 • Up to 2 or 6i^ values of x. may be handled in this Banner o 
For the particular value of x three-quarters of the way between 
,Tt « 3 X 2 and x.^. - .^ x 2 ^ x would "^.prenr in the AC rsj 

AC B-Register 



Jij ■- .:) 5C 



fflAx. = Oo75 X 2 



whe 



re,ix^ « 2** • Hlth this re? resentction of x In the AH, the Instruction 



gr *g (shift everythini^ In the AC anc B-Fef Ister 9 f IgH losition? to 
the right without roundoff) would leave x, x 2'^ In the AC, and 

m % r 'In the B-Eegletcr« Then x. 'x^'^^ nay be used to find f(x.) and 

-_ _-5^ ' a-Regl ster _ 
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^^^1*-!^* «0 was doM in Bectlon H«l* Th« si 15 Inetructlon (shift every- 
thi^ left 15 digit positions) would then put the fraction jjg in the AC* 
Thus the three values f(Xj), f(x.^.), and n, required to find f(x) from 
the extrapolation fomula, are all made available* 
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^Linear Interpolation in a Table of n Values (n^64) 





Tnp^rttrtji9fl 
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ad 203 
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4« 


td 11 


to find 


V* 


td 13 


fM. 




6. 


ad 200 


Jk 




7. 


td 10 


f(Xl+l),»* 




8« 


0l 15 




.9. 


ta 202 


rind f (x) 


/lO. 


^ - - - 




11. 


au ^ ^-^ 


by interpo- 


<12. 


nr 202 


lation 


13« 


ad 


fonsttla 


\M' 


ta 202 




15- 
1 


ep next Job 




t 






200* 


ri 1 




201. 


X 




202. 







203. 


ri 204 




204* 


f (x^(x 2*^^ 




205. 


f (x^ • 1) 




206. 
t 

■ 


f (x^ « 2) 




204' 


1 

fn. 


f (x^ « n) 



mw% 



X 

x^(x 2*^^) in AC, m x 2*^^ in B-Regiater 



addreae of f (x^) 



addreaa of ^(Xj^n) 
m in AC 



f(x,,,) 

f(x^^j) - f(xp 
[f (x^^^) . f (xj] 



f(x) « f (Xj) ♦ » [f(x^^j) ♦ f (Xj)] 



value of X for nhioh f (x) is deaireo 
receiTee n, then f (x) 
addreae of beginning of table 



tabulated fUnetiona 
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BStJl 



SiglM^ 



Approire^i 



AttaehMli A-36229 


la 


A-36230 


6a 


A-35676 


8a 


A-36237 


1^ 


A-3623it 


15a 


A-36232 


15b 


A-36233 


15o 


A-36231 


15d 
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The deeiocil Eyetra of nunberE is baied on the uee of 10 sopa* 
rate digits for oountlng* TboBo are the lotegere through 9* Corree- 
pondingly, the binexy syeten of onabere le baaed on the uee of 2 digite, 
axid 1* 

The muBber $097*23 in the decimaX iyaten nay be repreaented in 
Bore eaqplicit fonn, namely i 

5(10^) + O(IO^) * 9(10^) + 7(10°) ^ 2(10"^) + SdO*"^) 
The tuo forna are identioal except that each power of ten given in the latter 
fom ie ioplied in the f irat fonn by the number of plaoea to the left or 
right of the decinal point of ita ooeff ioient* Uae of the second fora of 
the nuBber ia hardly neoeasaxy in the decinal ayatev, but it ie helpful in 
explaining the binary ayatem, or aqy other nunber eysten* 
The following ia an exanple of a binary number t 
1110.101 
It alao may be repreaented in a eorreaponding escplicit form, uaing powera 
of 2 inatead of 10, and 2 poaeible ooeffloienta instead of 10, thuat 

1(2^) ^ 1(2^) + 1(2^) + 0(2°) ^ 1(2*^) +0(2*^) ♦ 1(2^^). 
Thia ie equivalent tos 

1(e) ♦ lU) + 1(2) + 0(1) * 1(1/2) + 0(1/A) + 1(1/8). 
Adding these factors together givea 14*625 aa the decimal representation 
of the original binary number, 1110* 101* 

Binary numbers in ^WI registers are either poaitive integers 
in the case of addresses in instructions, or less than one in magnitude 
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In th^ (Saee ®f nvm®rlaal ¥&lues to b® umd In calsulatione® In other 
words 9 the binary pr^int is eonslder«d to b© ats (l) the eaKtrssae right 
of the 11 digit positions a¥ailabl® for addresses j (2) the esdl^reme left 
of the 15 digits available for magnitudes of words representing immbers® 

ThuSy 98 in the esffisuple of the inore esqplicit form of represent* 
ing binary numbers, when the digit positions for an address contain a 
1 they represei^t the values indicated: 



,^^ 



,8 



25 2*l23i22 2^ 



Eangfi of addresses held in one register 
to 20A7 in steps of 1 

Cornispondingly for words representing tmmberss 




J- 



f'k' 



,-5 



■|»^$^ 



=8 



2-9|"^Q~^-r^=i^.i5| 



Ran|8 of positive §M negative nsjmbers 
to 1-2^^5 in steps of 2°°2.5 



etepj 

Sca]e faetoring permits the representation of a praeti©ally 
unlimited ran|ij ©f numbers* 



EstMple of an address im. a 3?©gi8t©r follows j 
The Bddress 14.9 san h® represented ass 
2*^ ->• 2^ + 2^ + 2° (s 128 -^ 16 -^ i4 > 1 ^ 149) 
In binary form this is 10010101« 
In E register the address 149 would bes 




i|o !oii|o|iTo7r 



order 5€cr/*/v 



2* 2^ 2° 
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SscBinple of musbers 

The number •*' «>375 can be represented asi 
2**^ + 2"^ (^ »250 + .125 « o375) 
This ie repreeented in a register ass 

FioirriTglciclcicic^oiOicicni^ 



The negatiire niunber -^ 9 375 is represented Iqr changing all the 

£ tc i-e and I'S t^ e of its positive representation: 

liliLMljUiilllliiili^^ 



Beruinse r«f this don7en1iion| the noaber ^^^a be repreeen^d as 

-^ with 16 0*8 or as -0 with 16 l»s* 
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B« DECIMAL TO BIMikRY CONVERSION 

The proceeding exanplee of an addreae and a nunbar In a register 
indicate a netbod for decimal to binary convereiono The addreee HS 
free converted to binazy form by first finding the largest power of 2 (2 ) 
which would fit into 14>9; next was found the largest power of 2 (Tr) which 
would fit into 21 y the difference between U9 and 2\ etc* It was shown 
that the sum of these powers of 2 equaled the original number to be con'* 
verted* Civen these powers of 2, the binary number was written out by 
placing I'^s in the corresponding binary digit position to the left of the 
binary point* 

The number 0375 was converted by the same process* Since 0375 
is less than 1, only negative powers of 2 fit into it, and these are rep*^ 
resented to the right of the binary point* 

A mechanized technique for handling this conversion for integer 
numbers ist 



2 IMS * iccioioi* 


2 74 1 >| 


^. 


2 37 




2 |18 1 




2 |9^ 




Z \U 1 




2 Li 




2 tl 




1 





Report R-»196 A v 

Coneidering thm binary point to be aboTo tbo first renAinder^ the result 
is 10010101 « H9» In this process the operator divides 2 into U9ff puts 
the integer result below, and the resainder as 1 to the right of fA* 
Next he divides 2 into lA, puts the inteier result below and sinoe there 
is no renainder puts($ to the right of 37* The process is continued until 
the last quotient is a 0« 

The sane thing is done with fraotions, except that the deciaal 
nunber and its resainder are divided bj 1/2 (sniltiplied by 2) each tiBe, 
and a 1 is used to indicate that the result of a step has exceeded 1, a 
to indicate the result of a step is less than le 

•one « jj2Z5 (x2) 

C(.750) (x2) 
1(.500) (x2) 
l(.CCO) (x2) 
O(.CCO) 

This is equal to eOllO, again considering the binary point to 
be above the first or lo 

To convert a binary to a decimal nuaber, sinply add the deciieal 
values of the powers of 2 represented by eaeh 1 in the binazy nuabero 
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fiKi bae ordere vhieh eblft the potition of a number either to 
the left or to the right in a register* 

A shift of one digit position to the left oorresponds to miti« 
plying the original number by 2* Shifting to the right corresponds to 
division by 2* This is so because the shift is made relative to a fixed 
binazy point* Correspondinglyp shifting digits in the decimal t^TStem 
relative to the decimal point corresponds to multiplication or division 
\^ 10* 



Report R496 A vli 

f Rpci^ yiflNG IN OCTAL (BAS E 8> FORM 

Binary=to«'Octal cottverelon is simpler and mor© rapid than la 
blnary"to«daoimal comrarelon* In the octal ayeten of nuinberlng,the eight 
dlglta through 7 are uaedo Three binary digits are sufficient to represent 
any single octal digit* 

000 

1 001 

2 010 

3 Oil 

4 100 

5 101 

6 110 

7 111 

The binary form of an octal rasmber may be found directly from the above 
table, which is short enough to be readily memorised* For exaaqple, the 
octal number 316 is Oil 001 llOo Any binary number »ay be concerted to 
octal by slowly reading off from the left the octal digits (each consisting 
of three binary digits)* 

Programs for WVI are often prepared in octal form rather than 
in decimal f orn because of this greater ease of converslono 

The VWX control panel has indicator lights which gi've the binary 

contents of any desired register* The binary register contents may be 

readily recorded in octal form or readily checked agsinst a program origin^ 

ally written in octal form* 

For programs written in octal form, registers are numbered 

C, 1, 2, 3, I, 5, 6, 7, 10, 11, 12, 13, U, 15, 16, 17, 20, 21, 

75 1 76, 77, !^00, 101* •» Special program forms numbered ootally are 
available at project Whirlwind* 
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COMPUTER PROGRAMS 

Program . A program is a sequence of actions by which a computer 
handles a problem. The process of determining the sequence of actions 
is known as programming. 

Flow diagrams . A flow diagram is a series of statements of what the 
computer has to do at various stages in a program. Lines of flow indi- 
cate how the computer passes from one stage of the program to another. 

Coded program . Programs andf low diagrams are largely Independent 
of computer characteristics, but instructions for a computer must be ex- 
pressed in terms of a code. A set of instructions that will enable a com- 
puter to execute a program is called a coded program, and the process 
of preparing a coded program is known as coding. 

Orders and operations . Individual coded instructions are known as 
orders and call'for specific operations such as multiply, add, shift, etc. 

The computer code . The computer code described here is that of 
Whirlwind I, an experimental computer using binary digits, single- 
address order code, parallel operation, and electrostatic storage. It is ex- 
pected that computers of this type will ultimately achieve an average 
speed of 50,000 operations per second. 

COMPUTER COMPONENTS 

Registers and words . A register has 16 digit positions each able to 
store a one or a zero. A word is a set of 16 digits that may be stored in a 
register. A word can represent an order or a number. 

Arithmetic element . Arithmetic operations take place in the arith- 
metic element, whose main components are three flip-flop registers, the 
A-register, the accumulator, and the B-register (AR, AC, BR). The 16 
digitpositionsof AR starting from the left are denoted by ARO, AR 1, . . ., 
AR 15. Similarly for AC, BR. Words enter AC through AR; BR is an 
extension of AC. 

Storage . The term "register" by itself refers to the main electrostatic 
storage, which consists of 2''-''- or 2048 registers, each of which is identi- 
fied by an address. These addresses are 11 -digit binary numbers from 
to 2047. The computer identifies a register by its address. 

Input-output. All information entering or leaving the computer is tem - 
porarily stored in the input-output register (lOR). The computer regu- 
lates the flow of information between the internal storage and lOR, and 
also calls for any necessary manipulation of external units. The descrip- 
tive names of the input-output orders were chosen for photographic film 
reader-recorder units, but the ordersare applicable to other types of ex- 
ternal equipment. 

Control element . The control element controls the sequence of com- 
puter operations and their execution. Instructions are obtainedfrom stor- 
age in the form of individual orders, each of which is represented by a 
single word. 

Inter-connections . The four main elements (storage, control, arith- 
metic, and input-output) are connected by a parallel communications sys- 
tem, known as the bus. 

REPRESENTATION OF ORDERS 

Operation section . When a word is used to represent an order the first 
(left-hand) 5 digits, or operation section, specify a nartlcular operation 
in accordance with the order code. 

Address section . The remaining 11 digits, or address section, are 
interpreted as a number with the binary point at the right-hand end. In the 
majority of orders this number is the address of the register whose con- 
tents will be used in the operation. In orders si, sr, the number specifies 
the extentof ashlft; inrf, rb, the number specifies an external unit; in rl, 
rs, the address section is not used. 

Example . The order cax has the effect of clearing AC (making all the 
digits zero) and then putting into AC the word that is in the register whose 
address is x. If q is a quantity in some register, the order needed to put q 
in AC is not caq but cax, where x is the address of the register that con- 
tains q. 



REPRESENTATIC») OF NUMBERS 

Single-word representations . When a word is used to represent a 
number the fir St digit indicates the sign and the remaining 15 are nimieri- 
cal digits. Forapositivenumberthesigndigitiszero,andthel5numeri- 
cal digits with a binary point at their left specify the magnitude of the num- 
ber. The negative -y of a positive number y is represented by comple- 
menting all the digits. Including the sign digit, that would represent y . (The 
complement Is formed by replacing every zero by a one and every one by 
a zero.) In this way a word can represent any multiple of 2"^5 from 
2~^^ - 1 to 1 - 2"-'-5 . Neither +1 nor -1 can be represented by a single 
word. Zero has two representations, either 16 zeros or 16 ones, which 
are called -t-O and -0 respectively. 

Overflow — increase of range and accuracy . With single -word 
representation the range is limited to numbers between 2-^5 - i and 
1 - 2"^' . Programs must be so planned that arithmetic operations will not 
cause anoverflowbeyondthisrange. The range may be extended by using 
a scale factor, which must be separately stored. Accuracy can be in- 
creased by using two words to represent a 30-digit number. 

COMPUTER PROCEDURE 

Sequence of operations . After the execution of an order the program 
counter in the control element holds the address of the register from which 
the next order is to be taken. Control calls for this order and carries out 
the specified operation. If the order is not sp or cp(-) the address in the 
program counter then increases by one so that the next order is taken from 
the next consecutive register. The sp and cp(-) orders permit a change 
in this sequential procedure. 

Transfers . A transfer of a digit from one digit position to another af - 
fects only the latter digit position, whose previous content is lost. 

Negative zero . The subtraction of equal numbers produces a negative 
zero in AC, except when AC contains +0, and -0 is subtracted from it. 

Manipulation of orders. Words representing orders may be handled in 
the arithmetic element as numbers. 

Procedure in the arithmetic element . The execution of an addition in- 
cludes the process of adding in carries; this process treats all 16 digits 
as if they were numerical digits, a car ry from AC being added into AC 15 . 
A subtraction is executed by adding the complement. Multiplication, divi- 
sion, shifting and round-off are all executed with positive numbers, com- 
plementingbeingperformed before and after the process when necessary. 
For round-off the digit tn BRO is added into AC 15. 

NOTATION FOR CODING 

Addresses . A coded program requires certain registers to be used for 
specified purposes. Theaddressesof theseregistersmustbe chosen be- 
fore the program can be put into a computer, but for study purposes this 
final choice isunnecessary, and the addresses can be indicated by a sys- 
tem of symbols or index numbers. 

Writing a coded program . Registers from which control obtains or- 
ders may be called action registers, and should be listed separately from 
registers containing other information, which may be called data regis- 
ters. A coded program iswritten out in two columns; the first contains 
the index number of each action or dat^ register, and the second column 
indicates the word that is initially stored in that register. In many cases 
part or all of a word may be immaterial because the contents of the regis- 
ter in question will be changed during the course of the program. This 
state cf affairs is indicated by two dashes, for example, ca— . 

The abbreviations RC, CR . Abbreviations used in referring to the 
register that contains a certain word or to the word in a certain register 
are 



RC 
CR 



. = (Address of) Register Containing . . . 
. = Contents of Register (whose address is) . 



The symbol ri x . When an address forms part of an order it is repre- 
sented by the last 11 digits of a word whose first 5 digits specify an opera- 
tion. An address x that is not t>art of an order is repre >Qiited by the last 
11 digits of a word whose first 5 digits are zero, whicii l-i ■ quivalent to 
specifying the operation ri. Thus the word for an unattiushed address. x 
may be written rl x. It could also be written x x 2"^' . 



